Hidden Ferronematic Order in Underdoped Cuprates 
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We study a model for low doped cuprates where holes aggregate into oriented stripe segments 
which have a vortex and an antivortex fixed to the extremes. We argue that due to the interaction 
between segments a state with macroscopic polarization is stabilized, which we call a ferronematic. 
This state can be characterized as a charge nematic which, due to the net polarization, breaks 
inversion symmetry and also exhibits an incommensurate spin modulation. Our calculation can 
reproduce the doping dependent spin structure factor of lanthanum cuprates in excellent agreement 
with experiment and allows to rationalize experiments in which the incommensurability has an order 
parameter-like temperature dependence. 

PACS numbers: 74.25.Ha, 71.28.+d, 75.25.-j 



The question whether there is a broken symmetry 
hidden order in the pseudogap phase of cuprate high- 
temperature superconductors is still a matter of de- 
bate. In this regard, a popular proposal is unidirec- 
tional charge [l|, or spin and charge order d i] known 
as stripes. 

Bulk evidence for long-range charge stripe order has 
only been seen in few compounds. |5|-|7[ La2-2:Sr2^Cu04 
(LSCO) shows low energy incommensurate spin scatter- 
ing, usually associated with some kind of stripe corre- 
lation, which undergoes a rotation towards the diago- 
nal direction below hole concentrations x < 0.05, where 
it even becomes elastic. (8l-[lo| The incommensurability 
d (proportional to the inverse of the spin modulation) 
evolves linearly S x from x = 0.02 up to the metallic 
phase where it saturates above x > 1/8 [ll| as expected 
for stripes. 
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Quasistatic incommensurate spin scattering along the 
Cu-0 bond direction has also been found in detwinned 
YBCO. 14, [3] In this case, the absence of signatures 
of long-range charge order and the evidence of rotational 
symmetry breaking l3 - lij points towards a nematic or- 
der. It has been proposed that this order arises either 
from incipient unidirectional fluctuating stripes [l3| or 
from an independent d-wave type nematic actor which 
preserves translational symmetry. |18| Other proposals 



for a broken symmetry state include orbital currents 
which break time reversal symmetry [Toj. spirals [l^ and 
a d-density wave. [2l|,[23 

In this paper we show that neutron scattering exper- 
iments in strongly underdoped cuprates can be under- 
stood in terms of a phase which breaks rotational and 
inversion symmetry. It is formed by oriented stripe seg- 
ments which do not need to have positional order thus 
we call it ferronematic. The segments are oriented be- 
cause they sustain a vortex and an antivortex of the an- 
tiferromagnetic (AF) order in the extremes (Fig. [Ij. Al- 
though formally the phase has zero ordering wave vector 
in the charge sector, we show that it induces magnetic 
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FIG. 1: (color online), a) Stripe segment for 8 holes (cor- 
responding to 4 VA pairs) obtained by minimizing the GA 
energy of on a 16 x 16 lattice in the one band Hubbard mod- 
els with t' /t = —0.2. The radius of the circles are proportional 
to the added hole density while the arrows are the staggered 
magnetization, b) Spin currents defined from the conservation 
of the z component of the magnetization [2^ . In the contin- 
uum limit the spin current is proportional to the gradient of 
the phase of the staggered magnetization. 



incommensurate peaks in excellent agreement with ex- 
periment (Fig. 2]). Remarkably, the order parameter is 
proportional to the incommensurability as suggested by 
experiment. [3, [3l 

We focus on non-superconducting underdoped LSCO 
which can be grown in a structure with only two twin 
domains with different population. Therefore similarly 
to the case of YBCO the one-dimensionality of the 
low energy (diagonal) spin response can be clearly re- 
solved, [allfl] 

for the (low) 



Previous variational computations |23l-l26 
doping induced spin structure in cuprates suggest the 
formation of magnetic vortex-antivortex (VA) pairs the 
relevance of which has been discussed e.g. in context of 
the destruction of long-range Neel order. In order to 
study these textures we have performed variational cal- 
culations based on the Gutzwiller approximation (GA) 
of the extended one-band Hubbard model. The ratio be- 
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tween on-site repulsion U and nearest-neighbor hopping t 
is set to U/t = 8 as suggested by previous studies. |28l.l29j 
The ratio between next-nearest neigbor hopping t' and t 
is taken as a material parameter. [30| 

Within this framework the ener geti cally most stable 
solution for two holes is a VA pair. [2J| Our GA compu- 
tations reveal that the interaction between the pairs is 
dominated by an anisotropic short-range core-core con- 
tribution which originates from the distribution of the 
localized holes. This results in a head-to-tail aggrega- 
tion of VA pairs which tend to form chains in such a way 
that only the vortex and the antivortex on the extremes 
contribute to the long range distortion. Fig. [Ija) shows 
the spin and charge structure for 8 holes corresponding 
to 4 VA pairs. Examination of the spin current (b) al- 
lows to visualize the VA pair nucleated at the extremes 
and the fact that the texture breaks inversion symmetry. 
Notice that the segments tend to form an antiphase do- 
main wall of the AF order although the transition from 
finite segments to infinite stripes is non trivial and will 
be discussed elsewhere. 
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Fig. [2] shows that segments formed by 2A^ holes have 
systematically lower energy than 2N spin polarons with 
the preferred orientation depending on parameters as dis- 
cussed below. Seen the segment as formed by N VA pairs 
we obtain the binding energy Ebmd between the pairs 
from E{N) — N *Ei + {N — l)Ebind- Parameters appro- 
priate for lanthanum cuprates (t' /t ^ —0.15 . . .— 0.2) [30|, 
yield a slight preference for diagonally oriented segments. 
In addition we have checked that in the low temperature 
orthorhombic (LTO) phase an anisotropy of t' /t along 
orthorhombic a- and b-axis favors the orientation of the 
VA-pairs along the a-axis. 

Until now we have neglected the long-range part of the 
Coulomb interaction which limits the infinite aggregation 
of VA pairs. In the spirit of Ref. [3l[ we may estimate 
its effect by considering the segments (each hosting Nc 
charges) embedded in a homogeneously charged back- 
ground. For large Nc the Coulomb energy per charge 
increases logarithmically with the number of holes. At 
a fixed doping concentration x the energy per planar Cu 
reads, 



EcHNc) + -/ + ^\EMnd\ 
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with Ec — / (eaao) a charging energy expressed in terms 
of the orthorhombic lattice constant Co, the static dielec- 
tric constant eo, and 7 = (Ei + Ebind)/2. The last term 
comprises the fact that shorter segments have less bind- 
ing energy. The energy is minimized by Nc = \Eund\lEc 
which leads to short segments of only few lattice con- 
stants at infinitesimal doping. Considering a Wigner 
crystal of segments and screening we find that the length 
of the segments tends to grow rapidly with doping and 
may even diverge at a critical doping, [s^ 
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FIG. 2: (color online). Energy of a chain of 2N holes cor- 
responding to TV V-A pairs aligned along the vertical (open 
symbols) and diagonal direction (full symbols) as compared 
to the energy of 2N polarons. The slope of the dashed lines 
corresponds to the binding energy between pairs. Parame- 
ters: U/t = 8, t' /t = —0.2. The upper right inset shows the 
binding energy as a function of t' /t for vertical (squares) and 
diagonal (circles) directions. Lines are guide to eyes. Com- 
putations where done in systems with up to 20 x 20 sites. 



It is useful to remember that a XY vortex can be 



mapped to a 2D Coulomb charge [37j or alternatively to 
current wires perpendicular to the 2D plane. The sign of 
the charges or the direction of the current is then deter- 
mined by the vorticity. Thus the VA pairs map into 2D 
electric dipoles in the Coulomb analogy or 2D magnetic 
dipoles perpendicular to the former in the current anal- 
ogy. In the latter case the magnetic field maps into the 
spin current which will be used below. 

We assume that quenched disorder will yield a state 
in which segments have no positional order. We have 
checked numerically [s^] that the combination of long- 
range dipolar interaction plus short range interaction fa- 
vors a ferronematic alignment of the dipoles. 

Information on the spin response of a large aggregate 
of segments from the GA is hampered by the finite size 
of the clusters. Since the textures are planar and we are 
interested in the large scale behavior we consider a clas- 
sical AF XF-model with nearest neighbor interaction J. 
The segments are modeled as a chain of vacancies which 
alternately correspond to the center of a vortex and an- 
tivortex. The spin structure is then determined from the 
minimization of the classical energy. In order to prevent 
annihilation of VA pairs we introduce a frustrating AF 
second next nearest neighbor coupling J' across the seg- 
ments. The value of J' can be obtained by comparing 
the phase change across a single segment with the GA 
calculation (Fig. [T]) which yields J' « J. 

As an additional parameter we have to fix the filling 
factor V which is the number of charges per segment 
length. Whereas the site-centered VA chain shown in 
Fig. [T] has J/ sa 1 there exists a similar plaquettc cen- 
tered structure with comparable (though slightly higher) 
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FIG. 3: (color online). Spin phase distribution for (a) macro- 
scopically and (b) randomly polarized distribution of stripe 
dipoles in a 160 x 160 sites system. Segments have a length of 
8 sites on the diagonal (4 VA pairs) and are represented by 
arrows indicating the length and polarization. With a filling 
factor V — 0.7 the total number of segment sites corresponds 
to a doping of x = 0.03. 



energy, which has i' sa 0.65. The filling factor of these 
structures is thus similar to those of (infinitely extended) 
diagonal site- and bond-centered stripes. [3J] The fol- 
lowing calculations assume v = 0.7. Note, however, that 
experimental data im ply an increase towards v = I close 
to the AF boundary. [10[ 

Fig. El^a) reports the spin phase distribution for a par- 
ticular distribution of stripe dipoles with segment length 
of 8 sites at a; = 0.03, all polarized along the [—1,1]- 
direction. One observes a monotonous increase on the 
phase of the staggered magnetization along the [1,1]- 
direction, although the distribution of stripe segments 
is completely random. In Fig. [3][b) , for the same distri- 
bution of segments, the associated dipole orientation is 
now completely random. In this case the system disag- 
gregates into large areas with similar phase. 

The numerical results can be understood by consid- 
ering the long-range distortion produced by the seg- 
ments. Using the magnetic analogy mentioned above 
and standard arguments from the theory of macroscopic 
dielectrics and magnetized systems [s^ it is easy to see 
that the system develops a macroscopic "magnetic field" 
equivalent to an average net spin current perpendicular 
to the segments orientation and proportional to the av- 
erage polarization. 

Analogously we can obtain this result by considering 
the influence of a collection of equally oriented segments 
on spins at point r 



S^r) = 5'oexp(iQ •r)cos$(r) 
S'«(r) S'oexp(iQ •r)sin$(r) 



(2) 



where Q = (tt, tt) is the AF wave- vector and the phase 
$(r) can be expressed by mapping the Cartesian plane 
into the complex plane {x,y) ^ z — x + iy, 



Eqs. (HI), ((3)) yield the angular spin distortion due to the 
Nseg segments with ± "charges" at Zi^± (vi^±). 

By direct integration we obtain that the average spin 
current in a square system of dimension LP' in the direc- 
tion perpendicular to the "dipole moment" Pi = r_|_ — r_ 
is (V<I>(r)) = Trpi x z/Lp with z pointing out of the plane. 
The macroscopic spin current is given by V3>(r)„iac = 
TTz X P where P is the macroscopic polarization in the 
charge analogy (P = pNseg/L'^ for identical dipoles). 
Clearly P plays the role of the ferronematic order pa- 
rameter. From Eq. ^ one obtains that the macroscopic 
spin current implies an incommensurate spin response 
perpendicular to the segments with q = V^{r)mac- This 
is one of our central results and relates the incommensu- 
rability with the ferronematic order parameter. 

The doping dependence of the order parameter can 
be estimated by considering segments (number of sites 
= Z + 1) which host Nc charges so that the filling factor 
ly = Nc/ (^-1-1). Since x = NcNgeg/Lp is the concentration 
of charge carriers one thus finds a linear dependence of 
the incommensurability on doping 
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$(z) = Im J2 - ^«:-) - - ^*,+)] ■ (3) 



It is interesting to remark that if one computes the in- 
commensurability of infinite stripes consisting of collinear 
domain walls (thus without spin current) with filling 
[1, [l2j it coincides with the incommensurability of fi- 
nite ferronematic segments with the same filling and at 
the same doping, i.e. such that the total length of the 
segments coincides with the total length of the stripes. 
Thus experimentally the two phases are not easily dis- 
tinguished in the magnetic channel and the main differ- 
ence will arise in the charge channel with equally spaced 
stripes producing Bragg peaks in contrast to diffusive 
scattering in the case of segments. 

For a set of configurations of macroscopically polar- 
ized VA segments as shown in Fig. [SJa) we now evaluate 
the magnetic neutron cross section for different dopings 
which in Fig. 2] is compared with elastic neutron scatter- 
ing data from Ref. 0. The specific scattering geometry 
(Fig. 2b of Ref. i) which is composed of two twin do- 
mains with population 2 : 1 has been taken into account. 
This gives rise to the asymmetry of the spectra since 
Qaf of the B-twin does not coincide with Qaf of the 
A-twin. The incommensurate peak position Qc is deter- 
mined by Eq. |4] and is essentially independent from the 
segment length / (i.e. the dipole moment p). [s^ On the 
other hand the length infiuences the peak width as can 
be seen in the lowest panel of Fig. Qj. By decreasing the 
dipole moment p at fixed doping, the increasing number 
of segments decreases the fluctuations of the dipole po- 
larization. The strength of the incommensurate response 
is then favored with respect to the commensurate one. 
As can be seen from Fig. 3] one finds excellent agreement 
with the experimental data for segments with 4 VA pairs 
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FIG. 4: (color online). Fits of the spin structure factor 
(LSCO) at different dopings for A'^c = 8 segments as ex- 
plained in the text. For x — 0.05 we also show spectra for 
Nc = 4:,6 segments for comparison. Computations have been 
done on lattices with up to 160 x 160 sites and we average 
over 20 — 30 segment configurations where the experimental 
resolution (horizontal bar in the lower panel, cf. Fig. 4 of 
Ref. 0) has been taken into account by convoluting with a 
gaussian. Data by courtesy of S. Wakimoto. 



(i.e. 8 sites) whereas shorter segments underestimate the 
intensity at Qaf (cf. lower panel of Fig. 

As shown in the inset to Fig. [5] \ Eiiind\ decreases with 
increasing \t' /t\, the reduction is, however, more pro- 
nounced for diagonal VA pairs. This results in a crossing 
of the alignments as a function of N for —0.3 < t' /t < 
—0.2. In addition the lattice structure can influence the 
orientation. This explains the fact that a diagonal incom- 
mensurate spin response has only been observed in non- 
superconducting undcrdoped lanthanum cupratcs where 
\t' /t\ is small |30| and which display a LTO structure in 
the underdoped regime. At higher doping an increased 
fraction of local octahedral tilts with LTT symmetry [s^ 
may drive a reorientation of the stripe segments from di- 
agonal to the copper-oxygen direction at the insulator 
metal transition x ~ 0.05. In Bi2Sr2CaCu208+5 com- 
pounds where t'/t is large [s^l we expect alignment of 
the segments along the Cu-0 direction resulting in an in- 
equivalence of hole density on the corresponding x- and 
-y oxygen sites. It therefore would be interesting to in- 



vestigate whether this feature can account for the intra- 
unit cell nematicity observed by scanning tunneling mi- 
croscopy [1^ in these materials. 

Apart from the excellent agreement found for the spin 
structure factor other facts point to the correctness of 
our interpretation. We expect that thermal fluctuations 
will disorder the dipole orientation reducing the order pa- 
rameter until a thermodynamic phase transition occurs. 
It is not clear whether the high temperature state will 
be only nematic as in Fig. [Sfb) with a second transition 
at higher temperature to an isotropic state of fully dis- 
ordered dipoles or if the transition will be directly to the 
isotropic state. In both cases the incommensurability 
should display an order parameter behavior. In agree- 
ment with this expectation an ordered parameter like 
temperature dependence of the incommensurability has 
already been noticed in YBCO. 14, IE] 



Since the ferronematic state breaks inversion symme- 
try one expects on general grounds [i^ that it will lead 
to a real ferroelectric distortion, i.e. to become multi- 
ferroic, through e.g. the inverse Dzyaloshinskii-Moriya 
mechanism. [41| Indeed, a ferroelectric state has re- 
cently been detected in strongly underdoped lanthanum 
cuprates. [s^l Unfortunately, a small number of free car- 
riers will make the effect undetectable with capacitance 
measurements which may explain why it has been seen 
only at very low temperatures. An appealing possibil- 
ity would be to look for inversion symmetry breaking 
with second harmonic generation (SHG) which docs not 
require perfect insulating behavior, [i^l We expect that 
the SHG signal as a function of temperature tracks the 
incommensurability. 

The macroscopic constant spin current can be consid- 
ered as an average spiral behavior of the spins. In this 
sense our theory has some similarity with the proposal of 
Ref. [2^. However, in our case the spiral is a collective 
effect which is slaved by the ferronematic order of the 
stripe segments. 

Concluding, we have proposed a new phase for strongly 
underdoped cuprates which breaks C4 rotational and in- 
version symmetry. Our theory provides a consistent ex- 
planation for the elastic incommensurate response seen 
by magnetic neutron scattering experiments and recon- 
ciles it with the lacking signatures of charge order. It re- 
mains to be seen how the nematic segments act as seeds 
which lead to smectic correlations (stripes) in some of 
the cuprate materials. 
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